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Abstract 



We present an ab initio full-potential linearized augmented plane-wave 
(FLAPW) study of the structural and electronic properties of the two bulk 
unstable compounds FeSi (CsCl structure) and FeSi2 (CaF2 structure) which 
have recently been grown by molecular beam epitaxy on Si(lll). We obtain 
equilibrium bulk lattice constants of 2.72 A and 5.32 A for FeSi and FeSi2, 
respectively. The density of states (DOS) of FeSi agrees well with experiment, 
and shows metallic behavior. In agreement with a previous calculation the 
DOS of FeSi2 shows a large density of d-states at the Fermi level, explaining 
the instability of the bulk phase. The electron charge distributions reveal a 
small charge transfer from Si to Fe atomic spheres in both compounds. While 
in FeSi the Fe-Si bond is indeed partially ionic, we show that in FeSi2 the elec- 
tron distribution corresponds to a covalent charge accumulation in the Fe-Si 
bond region. The reversed order of d-bands in FeSi with respect to FeSi2 
is understood in terms of crystal field splitting and Fe-Fe nearest neighbor 
cM-interactions in the CsCl structure, and a strong Si p/Fe d bonding in the 
fluorite structure, respectively. 
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I. INTRODUCTION 



The current interest in transition-metal silicides stems from both fundamental and tech- 
nological issues. For applications in microelectronics and optoelectronics the main focus 
is on silicides which can be grown epitaxially on Si such as, e.g., NiSi2 and CoSi2 whose 
small lattice mismatch with Si allows for the synthesis of high quality semiconductor/metal 
heterostructuresQ. These systems open the way for the study of intriguing fundamental 
questions, such as Schottky barrier formation, interface structure, growth kinetics, and sta- 
bility of novel epitaxial phases. Whereas the CaF 2 structure of both CoSi 2 and NiSi 2 are 
stable bulk phases of these compounds, FeSi 2 does not crystallize in the ffuorite form. At 
this stoichiometry and below ~ 1000°C its stable bulk form isiH /3-FeSi 2 , which is semicon- 
ducting. This orthorhombic structure can be regarded as a distorted ffuorite structure, and 
its greater stability with respect to the fluorite phase can be understood in terms of the 
electronic structure of the latter: a high density of d states (<i-DOS) at the Fermi level Ep 
leads to a Jahn- Teller- like instability of fluoriteB FeSi 2 . In /3-FeSi 2 the d levels are crystal 
field split at E F , and the opening of a gap! {E„ w 0.8 eV) leads to semiconducting behavior. 
At composition 1:1 in the Fe-Si phase diagramQ we find the compound e-FeSi. This structure 
is actually referred to as the "FeSi" structure, and its Pearson symbol is cP8. 

Under the growth conditions of molecular beam epitaxy (MBE) on Si(lll) the Fe-Si 
phase diagram is dramatically different^!: at a Fe/Si composition ratio of 1:2 not only 
epitaxial /3-FeSi 2 is found, but also the bulk unstable fluorite phase, denoted 7-FeSi 2 in Ref. 
|j. At composition 1:1 besides e-FeSi the CsCl structure of FeSi grows pseudomorphically 
on Si(lll). Hence we deal with two novel materials completely absent in the bulk phase 
diagram. Before turning to the electronic and structural properties of these two structures, 
let us briefly discuss the physical mechanisms pertinent to the occurrence of novel epitaxial 
phases. (We shall use the term "epitaxial phase" only for novel phases which do not exist 
in bulk form, in order to emphasize their distinctiveness from stable bulk phases which may 
also be present in epitaxial phase diagrams.) 

We distinguish two classes of epitaxial phases: («) epitaxially stable phases, where the 
strain conditions imposed by the substrate select a species with a lower epitaxial forma- 
tion enthalpy than the strained bulk phase at that particular substrate lattice constant!!; 
(m) epitaxially metastable phases, where kinetic barriers inhibit the transformation to a 
(epitaxially) more stable phase. 

For 1:1 FeSi, we argue that up to a critical film thickness the CsCl structure of FeSi 
belongs to class (z), i.e., it is epitaxially stable on Si(lll). Experimentally!] it is found that 
low temperature MBE growth of FeSi on Si(lll) yields coherent films in the CsCl structure 
up to ~ 100 A, above which the creation of misfit dislocations leads to partial relaxation. 
Films thicker than 15 A transform into relaxed e-FeSi upon annealing at ~ 300°C. The 
experimentally determined lattice mismatches with Si(lll) are +2% for FeSi (CsCl) and 
—6.4% for e-FeSi, respectively. We can rationalize these experimental findings as follows 
(see Fig. [I]): Due to the large lattice mismatch, pseudomorphic e-FeSi will have a very high 
strain energy, and a small critical thickness hi, whereas the much smaller mismatch between 
FeSi (CsCl) and the Si substrate makes this material more stable at small thicknessR Since 
e-FeSi is the stable bulk phase, however, above a critical thickness H c where it is completely 
relaxed, it is lower in energy than the CsCl structure at any thickness. Therefore, at some 
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thickness /i 2 the two energy curves cross, and the CsCl structure ceases to be the epitaxially 
stable phase. From experiment we infer h 2 ^5 H c pa 15 A. Due to kinetic barriers (e.g., an 
energetically favorable interface structure), the CsCl phase is observed in metastable form 
even above h 2 , and beyond its critical thickness h$ ~ 100 A up to ~ 1000 A. 

At composition ratio 1:2 the situation is less clear. Many epitaxial relations for (3- 
FeSi2/Si(lll) with different lattice mismatches are reported0, and a disordered Fei-^n^Si 
phase of CsCl symmetry with an iron vacancy concentration x = 0.5 acts as a precursor 
to the ordered CaF 2 phasei'0. We will show, however, that the CaF 2 phase is even better 
lattice matched to Si than the CsCl phase, and hence a similar picture as sketched in Fig. 
[I] may apply also in this case. 

Whereas the existence of epitaxial fluorite FeSi 2 is not that surprising in view of its 
relation to the bulk stable /3-phase and the existence of CoSi2 and NiSi 2 in this crystal 
structure, the CsCl structure of FeSi seems to be rather "exotic" . It is interesting to note that 
the only transition metal silicides with composition MSi which occur in the CsCl structure 
are those with0 M = Ru, Os, and Rh. Ruthenium and Rhodium are the transition metals 
situated in the Pd row of the periodic table just below Fe and Co, respectively. These 
silicides normally order in the "FeSi" structure type, i.e., the structure of e-FeSi, and only 
under certain growth conditions do they form the CsCl structureS. For instance, in the case 
of RuSi it was observed that the CsCl structure formed "more easily in mixtures that had not 
yet reached equilibrium" 0. We have argued above that FeSi films are epitaxially metastable 
above ~ 15 A, and only grow under the out-of-equilibrium conditions of low temperature 
MBE. These considerations suggest that although the CsCl structure is extremely rare for 
mono- silicides, it is nevertheless most likely to be formed by certain group VIII transition 
metals, among which we find Fe and Co. Indeed, in Villars' three-dimensional structural 
stability diagramB, both compounds are found in regions of coexistence of the MnP-FeSi- 
CsCl structure types, but almost on the boundary with the region occupied by the CsCl 
structure. Very recently, epitaxial growth of CsCl-type CoSi on Si(lll) has indeed been 
accomplished by MBEE3. The existence of CsCl-type NiSi, however, seems less probable. 

The great versatility of Fe to form a number of epitaxial silicides on Si — among which 
are found metals (FeSi(CsCl),7-FeSi 2 ), semiconductors or semimetals (/?-FeSi 2 , e-FeSi), and 
even magnetic compounds (Fe 3 Si, possibly 7-FeSi 2 ) — offers a wide spectrum of potential 
applications. Theoretical studies of the electronic and structural properties of the Fe silicides 
are therefore of great interest. Semiconducting /3-FeSi 2 has been studied by the linear 
muffin-tin orbital (LMTO)i, and by the augmented spherical-waveB methods. In Ref. f§, 
the electronic structure of the fluorite phase 7-FeSi 2 is also given. The magnetic compound 
FeaSi (DO3 structure) has been studied theoretically by Kudrnovsky et a/.0 We have recently 
published a brief account of the electronic band structure of FeSi in the CsCl structure, 
together with experimental results!. 

In this work we concentrate on the properties of the bulk unstable "high symmetry" 
compounds FeSi (CsCl structure) and 7-FeSi 2 (CaF 2 structure). We determine their lattice 
constants in order to better understand their existence as epitaxial phases on Si(lll), and 
give a detailed account of their electronic band structures. We focus on the bonding between 
Fe and Si both in FeSi and 7-FeSi 2 and compare to NiSi 2 and CoSi 2 . Special attention is 
given to the covalent versus ionic nature of the Fe-Si bond, and the role of Si p - Fe d 
interaction. 
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The paper is organized as follows. In Section [TJ we present the details of our ah initio 
calculations and the resulting equilibrium structural parameters of the two hypothetical bulk 
phases; the Fe-Si bond is studied in Sec. [TTI[ in Sec. [TV] the electronic band structures are 
discussed; and finally the conclusions are given in Sec. |V|. 



II. SELF-CONSISTENT CALCULATIONS 

FeSi is assumed in the CsCl structure (space group Pm3m), i.e., in a simple cubic lattice 
with Fe at the origin and Si at the body center. Both types of atoms are coordinated with 8 
nearest neighbors (nn) of the other species. 7~FeSi2 crystallizes in the CaF2 structure (space 
group Fm3m), i.e., in a face-centered cubic lattice with Fe at the origin, and two Si atoms 
at positions ±(-^~), respectively. Each Si atom is tetrahedrally coordinated with 4 metal 
atoms, and each metal atom has 8 Si nn. 

We have calculated the electronic energy eigenvalues, total energies, and wave func- 
tions of the two crystals with the self-consistent full-potential linearized augmented plane- 
wave (FLAPW) method^, neglecting spin-polarization. The exchange-correlation energy is 
treated in the local density approximation (LDA), using the Hedin-Lundqvist prescription. 
Wave functions are expanded in spherical harmonics with I < 8 (I < 4 for total energy 
minimization) inside the atomic spheres whose radii are 1.85 a.u. and 2.3 a.u. for Si and Fe, 
respectively. In the interstitial region a plane wave basis is used with wave vectors up to 
3.6 Ry (8.5 Ry) for wave functions (charge and potential). Integration in fc-space is approx- 
imated by summing over up to 48 fc-points in the irreducible wedge of the corresponding 
Brillouin zone. The energy band structure and density of states (DOS) of the equilibrium 
structures are evaluated using higher cutoffs and more fc-points. 

The lattice constants and bulk moduli obtained by total energy calculations are summa- 
rized in Tab. |, together with the lattice constants of the bulk phases. We have fitted the 
total energies calculated at six different lattice constants to both Murnaghan's equation of 
stateEJ and a third- or fourth-order polynomial. The predicted equilibrium lattice constant 
is independent of the fitting procedure, whereas the fluctuations of the fitted bulk modulus 
provide the theoretical error bars given in the Table. The theoretical values predict a virtu- 
ally perfect lattice match of CsCl-type FeSi and Si. Recent X-ray diffraction datJ on thin 
epitaxial FeSi films yield a lattice constant a^gj = 2.77 A, which is 1.8 % larger than our 
theoretical value, and a larger lattice mismatch with Si, i.e., 2apgg ; = 5.54 A (a^ p = 5.43 
A) which remains however rather small. The predicted lattice constant of 7-FeSi2 is only 2 
% smaller than that of FeSi (Tab. |), even though there are only half as many iron atoms 
per mole as in FeSi. Assuming the same magnitue of error on our predicted lattice constant 
as for the CsCl structure, we obtain a corrected value of 5.42 A, which is almost perfectly 
lattice matched with Si. This explains why the fluorite phase can compete with the lattice 
mismatched bulk phase (/3-FeSi 2 ) in the presence of the Si substrate (see our discussion in 
the Introduction). While the Fe-Si nn distance differs by only 2 % in the two crystals, the 
Fe-Fe nn distance is 2.72 A in FeSi, and 3.76 A in FeSi2- Note that this Fe-Fe spacing in FeSi 
is only 10 % greater than the nn distance, and 5 % smaller than the next-nearest neighbor 
distance along the cube edge in bcc iron. 

Given the near lattice matching of both FeSi (CsCl) and FeSi 2 (CaF 2 ) with Si, we shall 
describe the electronic properties of the unstrained hypothetical equilibrium phases, rather 
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than the rhombohedrally distorted ones observed on Si(lll). The strain induced changes 
on the electronic band structure are expected to be very small. 



III. THE Fe-Si BOND 

The nature of the chemical bond in transition metal silicides has been the subject of 
extensive theoretical and experimental studies&H. In general, the bonding in the silicides 
is understood in terms of metal d - silicon p hybridization. There has been a controversy 
about the bonding in NiSi2- Tersoff and HamannS proposed that the bonding in NiSi2 and 
CoSi2 is similar to that in Si, an idea first put forward by Chabal et a/Jll Following their 
arguments, the tetrahedrally coordinated Si atoms build sp 3 hybrids, which in the silicides 
form directed covalent bonding states with Fe d orbitals. In a later study, Lambrecht et alB^ 
came to a similar conclusion, but attributed more importance to the charge transfer from 
Si to Ni which contributes an electrostatic term to the cohesive energy. It is well known, 
however, that the charge transfer is not uniquely determined, since it depends on the choice 
of the "atomic spheres" where the charge is measured. With this in mind, we compare in 
Tab. the valence charge contained in the atomic spheres used in our calculations. Besides 
the self-consistent density of the two crystals we have calculated the charge density of the 
isolated atoms as well as the lattice superpositions of the atomic densities. In FeSi 2 , the 8 
Si nn of Fe contribute more charge to the Fe spheres due to the smaller Fe-Si distance. On 
the other hand, the six Fe nn of Fe in FeSi are much closer than the 12 Fe nn of Fe in FeSi2, 
therefore superposition of atomic densities gives a Fe charge 0.08 electrons larger in FeSi than 
in FeSi2- As to the charge in the Si spheres it is obvious that twice the number of Fe nn of Si 
in FeSi yields a larger overlap than the corresponding one in FeSi2- The remaining free-atom 
charge is distributed in the interstitial region. When comparing the free-atom superpositions 
with the self-consistent densities, we observe an additional net charge transfer AQ-p e of 0.18 
electrons into the Fe sphere in FeSi, and one of 0.35 electrons into the Fe sphere in FeSi2. 
Note that now the Fe sphere in FeSi2 contains more charge than the one in FeSi. It would 
be misleading, however, to draw conclusions on the ionicity from these numbers alone, due 
to the ambiguity pointed out above. We will show that the excess charge in the FeSi 2 Fe 
sphere with respect to the FeSi one is due to a piling up of covalent charge in the Fe-Si bond 
region of FeSi2- 

To this end we analyze the charge densities instead of the integrated charge in the 
spheres. In Fig. |2| we display contour plots of the difference between the self-consistent 
charge density and the superposition of atomic densities in the (110) plane of FeSi2 (a) 
and FeSi (b). Besides the obvious difference due to the presence of an additional Fe atom 
in FeSi, we notice a weak piling up of charge midway between Fe and Si in the ffuorite 
structure. This covalent charge accumulation is almost entirely contained in the Fe sphere, 
and is responsible for the larger charge transfer in FeSi2 discussed above. In FeSi, on the 
other hand, the charge re-distribution is unimportant in the interstitial and bond region, 
and we interpret the charge transfer of 0.18 electrons from the Si to the Fe spheres as a real 
ionic character of the bond. This conclusion will be reinforced by the crystal field arguments 
presented in the next section. 
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IV. ELECTRONIC STRUCTURE 



The energy dispersion relations of 7-FeSi2 and CsCl-type FeSi can be found in Refs. |] 
and §, respectively, and are not reproduced here. In Fig. |3| we show the DOS of the two 
materials instead^. In Tab. |T| we provide the angular-momentum projected densities at 
Ep and integrated up to Ep, respectively. Our DOS for FeSi2 is in perfect agreement with 
the LMTO one by Christenseni. The two main peaks in the valence spectrum originate 
from the metal d/Si p bonding states at ~ —4 eV (r^, states at the zone center) and 
from the non-bonding metal <i-band manifold centered at ~ —1.5 eV (r^ 2 states at the 
zone center), respectively. LDA energy eigenvalues at high symmetry points are listed in 
Tab. |TV| , together with the angular-momentum decomposition of the wave functions inside 
the atomic spheres. The pd hybridization broadens the metal d bands considerably as can 
be seen from the DOS and the orbital characters of the eigenfunctions. Given the strong 
similarity between the present Si s and p partial DOS with those of CoSiJ^I, we conclude 
that the binding mechanism in FeSi2 is basically the same as in the former compound. It 
has been shown by Christensen that the rigid band assumption is nearly fulfilled for the 
three compounds FeSi2, CoSi2 and NiSi2 (see Fig. 1 in Ref. |]), and therefore the number 
of valence electrons pins the position of the Fermi level in an otherwise similar DOS. At 
the Fermi level there is a rather strong peak in the <i-DOS of FeSi2 (28.5 Ry -1 , cf. Tab. 
IIJ), with an appreciable admixture of Si p-states. The location in fc-space of these states is 
shown in Ref. [| to be on the hexagonal Brillouin-zone face. We have previously estimatedi 
that the Stoner factor is considerably larger than unity and explained bulk instability of 
the fluorite phase in terms of the Stoner model. Indeed, a spin-polarized calculation by 
Christensen predicts a ferromagnetic moment of 0.3 /ie per formula unitS. Experimental 
confirmation of ferromagnetism in 7-FeSi2 is still lacking, since the grown films are very thin 
(< 15 A) and magnetic measurements are difficult to perform. Above Ep we observe a Si 
sp antibonding peak at 3 eV, which seems to be even more pronounced than in CoSi 2 0, 
and which is a typical fingerprint of sp 3 hybridization. Together with the discussion of the 
charge transfer in Sec. |T| the picture of covalent bonding in fluorite FeSi2 is thus confirmed, 
similarly to the case of CoSi2 and NiSi2, and in agreement with the interpretation of Tersoff 
and Hamannlll for the latter two compounds. In FeSi2, however, the covalent contribution 
to the free energy is too small to stabilize the fluorite phase in bulk form. The kinetics 
and epitaxial constraints of iron silicide grown by MBE on Si(lll), however, lead to the 
formation of fluorite FeSi2 as an intermediate step of the disordered Feo.5no.5Si to /3-FeSi2 
transition!! at ~ 500°C. 

Let us now turn to FeSi (right panel in Fig. |3|, and Tab. [V]). Again the total DOS is 
dominated by the metal d bands, giving rise to a strong 2 eV broad peak centered at ~ —2 
eV. The Fermi level is located at the upper end of a region of low DOS, at the onset of 
a less pronounced, sharp peak of Fe d states. The valence DOS agrees very well with the 
ultraviolet photoemission spectrum of Ref. |5], provided that the theoretical curve is shifted 
by ~ 1 eV towards E F . This final-state effect is discussed in more detail in Ref. ^} A very 
similar DOS is also predicted for CsCl-type RuSi by a LMTO Green's function calculation of 
Ivanovskiiii. There it is also shown that RhSi has a virtually identical DOS, except for the 
position of Ep which now falls at higher energy in the sharp d peak. The same behavior can 
therefore be anticipated when replacing Fe by Co in FeSi, indicating that the hypothetical 
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CsCl phase of CoSi is even less stable than FeSi, and likely to be magnetic. Below Ep the Si 
s and p projected DOS is rather similar to that in FeSi2- This does not imply the formation 
of sp 3 hybrids, however, since Si is now 8-fold coordinated, instead of 4-fold as in the ffuorite 
structure. We attribute the relatively large width of the Si s band to bonding with Fe s and 
p states (see Tab. [Vp. 

An important difference between the band structures of the two materials lies in the order 
of the Fe d bands at T: in FeSi T25' lies above T^, whereas in FeSi2 T25' lies below T^. An 
inverted crystal field splitting was previously noted in band calculations of chalcopyritesil 
(CuInSe2) and in II- VI compoundsc3. The results were explained in terms of symmetry- 
enforced pd couplings. A similar interpretation is pertinent here: In FeSi the Si p-bands 
(x, y, z) transform like T^, and thus do not mix with the d-band triplet (xy, yz, zx) of Fe, 
which transforms like ^5'. Hence, the crystal field splitting corresponds to Ti2 below T25' like 
predicted by the point-ion model. This mixing is allowed, however, in the ffuorite structure, 
where the p-orbitals of the two Si atoms in the unit cell originate both T 15 and T 2 5' states. 
The resulting pd bonding lowers the T 2 5' level relative to the non-bonding T 12 level, thus 
reversing the sign with respect to that originating from the crystal field splitting alone. This 
effect is illustrated in Fig. [|, where the energy levels of the Fe d and Si p states are shown for 
FeSi2 in the CaF 2 structure, (a) without and (b) with taking the pd interaction into account. 
The pd repulsion found in certain II- VI zincblende semiconductors has a similar origin!!: the 
cation d bands of r 15 symmetry mix with the anion p bands and are pushed below the T 12 
d states, opposed to crystal field considerations alone. Note, however, that the pd repulsion 
in ZnTe, CdTe, and HgTe, e.g., is much weaker than in the present case of FeSi2 and in the 
case of the chalcopyrites, since (i) the cation d - anion p energy separation is larger for these 
atom pairs than for Fe and Si, {it) the cation-anion distance is considerably larger than the 
Fe-Si one in FeSi2- Therefore, the r 15 — T 12 splittings are much smaller in magnitude in 
these II- VI compounds!! than the T25' — Ti2 splitting in FeSi2, which we predict to be —2.55 
eV. 

In FeSi, on the other hand, E{T<2&) — E(T 12) = +3.24 eV, which is an unusually large 
crystal field splitting. We argue that it has two major contributions, an ionic one and 
one due to repulsion between Fe d states of neighboring Fe atoms. Placing positive point 
charges +ze on the Si atom sites and negative charges — ze on the Fe sites will lower the 
electrostatic binding energy of the Ti 2 states relative to the T 2 5' ones. We estimate the 
ionic contribution AEi to the splitting by assuming analytic Slater orbitals for the Fe 3d 
wave functions, and by performing the lattice sum involved in the fourth-order non-spherical 
portion of the crystal potential as described by Callaway!! for cubic crystals. We obtain 
AEi ~ +1.14,2 eV for Fe _z Si +z . Clearly, this crude point-charge ionic model is not sufficient 
to quantitatively explain the level splitting, but the important points are that the sign agrees 
with the calculated one and that ionicity contributes to the splitting. A similar calculation 
for the fluorite structure (placing the charge — 2ze on the Fe sites) yields a d level splitting 
with equal sign as in FeSi, but 40 % smaller. This is not surprising, since the nearest neighbor 
configuration of the Fe atoms is the same for both structures, and the contribution from the 
outer shells decreases rapidly^. Hence, neglecting p-d bonding in the fluorite structure the 
Ti2 states would be lower in energy than the T25' ones (Fig. 01(a)), very much like in the 
CsCl structure. This shows once more that ionicity is not important in FeSi2- 

The other effect which contributes to the large splitting in FeSi is the Fe-Fe interaction. 
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We noted in Sec. [D] that the Fe-Fe nn distance in FeSi is rather close to the equilibrium 
bond length in bcc iron. According to Harrison! the dd interatomic matrix elements decrease 
with the 5th power of the distance, and should therefore be about 60 % as large as in bcc 
iron. Recent tight-binding fits to the present band structures show indeed that the Fe-Fe 
dd interaction cannot be neglected in FeSi, whereas it is less important in! FeSi 2 . In FeSi 
the Fe atoms occupy simple cubic lattice points, and hence are octahedrally coordinated 
to their nearest Fe neighbors. The Ti2 orbitals form a-bonds, whereas the r 2 5> orbitals are 
7r-bonded. Using Harrison's universal tight-binding parametersil a splitting of ~ 2.2 eV 
is obtained. Again, this is not sufficient by itself to explain the splitting obtained in the 
self-consistent calculation. Considering in addition the ionic contribution discussed above, 
however, (calculated with a better approximation than the point-charge model) one should 
recover the r 25 /-ri 2 splitting of the FLAPW calculation. 

Inspection of Tab. |V| and of Fig. |3] (right panel) shows that Si p/Fe d hybridization is 
present also in FeSi, although to a smaller extent than in FeSi2, at fc-points other than T. 
At the .R-point, e.g., the fc-dependent phase factors associated with the atomic orbitals in 
a LCAO scheme cause the Si (x, y, z) - like and the Fe (xy, yz, zx) - like wave functions 
to belong to the same irreducible representation i?i 5 , as they do under the Ta point group. 
Also along the A-line pd hybridiztion is allowed by symmetry, and leads to anticrossing (see 
Fig. 4 in Ref. |5|) of the twofold degenerate A 3 bands originating from both T 12 and r 25 / 
(retaining only dd matrix elements, these bands would cross along A). Due to prf-bonding 
the i?i5 state is pushed below Ti 2 , and contributes to the peak in the DOS at ~ —6 eV (see 
Fig. H). The upper A 3 band ends in a "drf-antibonding" R\ 2 state close to the Fermi level, 
with 99 % c?-character in the Fe atomic spheres. Due to the strong dd interaction we refrain 
from using the term "non-bonding d" for these states. 

V. CONCLUSIONS 

We have presented the first theoretical study of the electronic and structural properties 
of the bulk unstable cubic iron silicide FeSi(CsCl), and have compared them with those of 
7-FeSi 2 (CaF 2 ). Both structures exist epitaxially on Si(lll) due to a lattice mismatch with 
Si wich is smaller than that of the bulk stable phases e-FeSi and /?-FeSi 2 , respectively. In 
addition, kinetic barriers are believed to exist which account for the observed metastability 
of FeSi (CsCl) films at thicknesses exceeding its epitaxial stability range (Fig. [I]). 

We have established the close resemblance of the bonding configuration of FeSi 2 with 
CoSi 2 and NiSi 2 , i.e., Si sp 3 hybridization, strong metal d/Si p bonding, the presence of 
directed covalent Fe-Si bonds, and negligible ionicity. The giant density of d states at the 
Fermi level leads to a Jahn- Teller instability of the bulk compound, explaining its absence 
in the bulk phase diagram. 

CsCl-type FeSi, on the other hand, is characterized by a small charge transfer from Si 
to Fe, strong Fe-Fe dd bonding, and to some extent also by Si p/metal d interaction. We 
have shown that both the charge distribution and the electronic band structure of the two 
materials are consistent with the bonding mechanisms presented in this work. An empirical 
tight-binding study of the Fe silicides is under way!, and will address the lattice dynamics 
and its implications on the 7 — ► /3 phase transition in epitaxial FeSi 2 /Si(lll). 
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Comparison with other known monosilicides in the CsCl structure has allowed us to spec- 
ulate on the possible existence of epitaxial CoSi(CsCl) on Si(lll). Subsequent experiments 
have confirmed the epitaxial stability of this phase. A rigid-band model implies a rather 
high d-DOS at the Fermi level, which may cause the compound to order magnetically. 

The disordered CsCl phase Fei-^D^Si is observed over a wide range of Fe vacancy con- 
centrations < x < 0.5. From the present study of the ordered end compounds FeSi and 
FeSi2 we presume that the ionic, or Madelung contribution to the cohesive energy (present 
in the CsCl, but not in the CaF 2 phase) stabilizes the epitaxial disordered Feo.5Do.5Si CsCl 
structure relative to the ordered CaF 2 structure. In the latter the Fe atoms are located on 
alternate {111} planes, and the Si atoms are tetrahedrally coordinated with the Fe atoms. 
The role of the d electrons, however, cannot be obtained from the properties of the end 
compounds alone. A rigid-band model would predict the Fermi level to fall in a region of 
high d-DOS for increasing x, which is hardly reconciled with the experimentally observed 
epitaxial stability and absence of magnetic order in the defect CsCl phase. Further the- 
oretical studies using, e.g., the coherent phase approximation, are needed for a complete 
understanding of the electronic properties of Fei_ x D x Si as a function of x. We hope that 
the present study, which adds to the general understanding of transition metal silicides, will 
stimulate further theoretical and experimental work on the fascinating iron silicide system. 
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FIGURES 



FIG. 1. Schematic diagram of eptiaxial energy versus film thickness for FeSi on Si(lll). The 
CsCl structure is more closely matched to the Si substrate than the bulk stable e-phase, and hence 
is epitaxially stable at small film thickness. Above its critical thickness H c e-FeSi is the most stable 
phase, and the occurrence of the CsCl phase in this range can be explained by the existence of a 
kinetic barrier. Upon annealing at ~ 300° C the CsCl phase transforms irreversibly into the bulk 
e-phase. 

FIG. 2. Contour plot of the difference between the self-consistent charge density and the su- 
perposition of atomic densities in the (110) plane of CaF2-type FeSi2 (a) and CsCl-type FeSi (b). 
Areas of negative values are shaded, contours are labeled in units of 10~ 2 electrons/ (bohr) 3 , and 
the level spacing is 0.02 in these units. AQ a is the difference between the self-consistent valence 
charge and the superposition of the atomic valence charges in the sphere centered on atom a. The 
Fe-Si bonds are indicated with solid lines. 

FIG. 3. DOS of CaF 2 -type FeSi 2 (left panel) and CsCl-type FeSi (right panel), (a) and (c) 
show the total (solid line) and partial Fe d (dashed-dotted line) DOS in eV _1 per unit cell, (b) and 
(d) show the partial Si p (solid line) and Si s (dashed line) DOS. The partial DOS are projected on 
one atomic sphere in both structures. The energy is measured relative to the Fermi level (indicated 
by a solid vertical line). 

FIG. 4. Energy levels of Fe d and Si p bands at T in the CaF 2 structure of FeSi 2 , (a) neglecting, 
(b) including Fe d/Si p interaction. Note that pd repulsion pushes the T 2 5/ states below the Ti 2 
states, whereas crystal field splitting alone (a) predicts the Ti 2 states to be lowest. 
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TABLES 

TABLE I. Structural definition of the stable bulk phases! of FeSi and FeSi2, as well as of the 
bulk unstable phases which exist epitaxially on Si(lll). For the latter the theoretical predictions 
of the present work are compared to experimental values, where available, and to other theoretical 
results. The experimental lattice constant of Si is 5.43 A. 



Compound 


Space group 




Lattice parameters (A) 


Bulk modulus (Mbar) 


/ PpQrcnn ^^rinhnl 1 






(exp.) 


(theo.) 


(thpn 


Rulk stable- 
















CI 


= 4.46 






(cP8) 












/3-FeSi 2 


Cmca 


a 


= 9.86 






(oC48) 




b 

c 


= 7.79 
= 7.83 






Bulk unstable: 












FeSi (CsCl) 


Pm3m 


2a 


= 5.54 a 


2a = 5.44 


2.70 ±0.07 


(cP2) 












Fei-^D^Si 


Pm3m 


2a 


= 5.40 a 






(CsCl, x < 0.5) 




(x 


= 0.5) 






7 -FeSi 2 (CaF 2 ) 


Fm3m 






a = 5.32 


1.90 ±0.05 


(cF12) 








a = 5.39 b 


2.06 b 



a Reference |^. 
b Reference 



TABLE II. Valence charge contained in the atomic spheres of radius R and in the interstitial 
volume of CsCl- type FeSi and CaF 2 -type FeSi 2 . The total valence charge of Fe is 8, and that of Si 
is 4 electrons. 





Charge 


Fe 

R = 2.3 a.u. 


Si 

R = 1.85 a.u. 


Interstitial 


Total 


Fe, Si 


isolated atoms 


6.24 


1.25 






FeSi 


free-atom superposition 


7.30 


1.66 


3.04 


12.0 




self-consistent 


7.48 


1.48 


3.04 


12.0 


FeSi 2 


free-atom superposition 


7.22 


1.55 


5.68 


16.0 




self-consistent 


7.57 


1.42 


5.59 


16.0 
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TABLE III. Partial number of states (NOS) below the Fermi energy Ef, and density of states at 
Ep (DOS in Ry _1 ) normalized to one atomic sphere (cf. Tab. ||) in CsCl-type FeSi and CaF2-type 



FeSi 2 . 







Fe s 


Fe p 


Fe d 


Si s 


Si p 


Si d 


FeSi 


DOS 


0.001 


0.018 


1.678 


0.051 


0.020 


0.004 




NOS 


0.385 


0.530 


6.528 


0.599 


0.720 


0.143 


FeSi 2 


DOS 


0.317 


3.740 


28.478 


0.073 


2.115 


0.449 




NOS 


0.466 


0.540 


6.544 


0.607 


0.706 


0.099 



TABLE IV. FLAPW energy eigenvalues at high symmetry points of the face centered cubic 
Brillouin zone for FeSi2 in the CaF 2 structure. Energies are measured relative to the Fermi level. 
The probability of the orbital character of the wavefuntions inside the Fe and Si atomic spheres is 
given in percent. In the last column we list the probability in the interstitial region. 



btate 


Energy 
(eV) 


s 


le 
P 


d 


s 


Si 

p 


d 


int. 


1 1 


-13.0 I 


1 A 

14 






on 
30 






00 


1 25' 


-3.95 






50 




17 




32 


■pi; 
1 12 


-1.40 






80 






3 


17 


-pc 
1 2' 


0.87 








67 






29 


pc 
1 15 


3.85 




28 






2 


13 


56 


x%, 


-8.93 




10 




36 




1 


52 


x\ 


-7.09 


15 




8 




19 




58 


xi 


-6.43 






28 


33 




2 


37 


xt 


-3.49 




13 






27 




59 


x\ 


-1.44 


14 




52 




1 


5 


28 


xi 


-1.43 






86 




4 


1 


9 


xi 


-0.67 






91 








8 


x\ 


2.60 






74 


16 






10 


L{ 


-9.73 


13 




5 


26 


5 




51 


T v 


-9.16 




10 




26 


6 




58 


LI 


-3.41 






53 




13 


1 


32 


L\ 


-1.70 


9 




63 


2 


5 


2 


17 


J v 

L 3 


-1.04 






87 






1 


11 


T v 


-0.24 




16 






26 


4 


53 


T c 
Li 2 , 


2.74 








23 


11 


3 


59 
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TABLE V. FLAPW energy eigenvalues at high symmetry points of the simple cubic Brillouin 
zone for FeSi in the CsCl structure. Energies are measured relative to the Fermi level. The 
probability of the orbital character of the wavefuntions inside the Fe and Si atomic spheres is given 
in percent. In the last column we list the probability in the interstitial region. 



State 


Energy 




Fe 






Si 




int. 




(eV) 


s 


p 


d 


s 


P 


d 




T v 


-14.60 


24 






27 






48 


-pv 
12 


-4.12 






76 






3 


20 


-pv 

1 9c/ 


-0.88 






95 






1 


3 


r c 

15 


3.38 




17 






37 




40 




-6.30 






50 




16 




34 


r>v 
riy 


-1.27 








64 






30 


R 12 


0.07 






99 








1 


R 15 


1.67 




39 








12 


48 


YV 


-10.47 




17 




34 






48 


x\ 


-9.04 


26 




10 




17 




46 


xi 


-3.15 






87 








12 


xi 


-2.53 






81 




6 




11 


x\ 


-0.69 






97 








2 




0.34 


4 




89 
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